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(57) Abstract: A method is disclosed for generating images of a human or animal 
body on the basis of 3D-reconstructions from 3D-XRAY or 3D-Computer Tomogra- 
phy measurements, which bodies comprise both natural tissue and one or more high- 
density objects. The method comprises the steps of executing the measurements, dis- 
tinguishing the one or more high-density objects and executing a separating procedure 
thereon for generating an improved image of regions of the natural tissue. In partic- 
ular, the method comprises the following steps: applying a ramp filter on the various 
projection measurements to single out the one or more high-density objects;segment- 
ing the singled-out one or more high-density objects into a separate 3D reconstruc- 
tion; suppressing the reconstructed one or more high-density objects from the original 
projection measurements ;and segmenting said projection measurements without the 
suppressed one or more high-density objects. 



wo ^7017263 A2 llllillllllillllllillil([illlllilllllil 



(84) Designated States (regional): ARIPO patent (GH, GM, 
KE, LS, MW, MZ, SD, SL, SZ, TZ, UG, ZM, ZW), 
Eurasian patent (AM, AZ, BY. KG, KZ, MD, RU. TJ, TM). 
European patent (AT, BE, BG. CH. CY, CZ, DE, DK, EE, 
ES, FI, FR. GB, GR, HU, IE, IT, LU, MC, NL, PT, RO, 
SE, SI, SK. TR), OAPI patent (BF, BJ, CF. CG, CI, CM. 
GA, GN, GQ. GW. ML, MR, NE. SN. TD, TG). 



Published: 

— without international search report and to be republished 
upon receipt of that report 

For two-letter codes and other abbreviations, refer to the "Guid- 
ance Notes on Codes and Abbreviations" appearing at the begin- 
ning of each regular issue of the PCT Gazette. 



wo 2004/017263 




T/IB2003/003351 



A method and apparatus for generating an improved image of natural tissue in reconstructing 
body images from SD-measurements 



BACKGROUND OF THE INVENTION 

The invention relates to an apparatus for generating images of human or 
animal bodies on the basis of 3D-reconstructions from 3D-XRAY or 3D-Computer 
Tomography measurements, which bodies comprise both natural tissue and one or more 

5 high-density objects, such as being ftirthermore recited in the preamble of Claim 1 . A 
principal category of such objects is represented by intentionally introduced objects for 
maintaining or improving the quality of human or animal life, such as objects being in the 
form of surgical implants made from metal or other substances, reconstruction screws, plugs 
filled into teeth, coils introduced into blood vessels, and various others. A secondary category 

10 is without limitation formed by high-density markers used for allowing a registration to 

match various different data sets. The 3D reconstructing methods recited supra, and possibly 
others as well, are suffering from the visual artifacts that such high density objects may cause 
in their neighbourhood, and which artifacts will lessen the quality of the eventual image, and 
thereby diminish its value for diagnostic, curative and other purposes. In consequence, it 

1 5 would be advantageous to have an approach for suppressing such artifacts. The inventors 
have recognized the advantage of suppressing the high-density object(s) from the processing 
in an early stage of the latter. 

SUMMARY TO THE INVENTION 
20 In consequence, amongst other things, it is an object of the present invention 

to provide an image generating method that does not suffer from the above artifacts from the 
high-density bodies introduced. 

A further object of the present invention is to supplement the image from the 
. natural tissue in a secondary processing stage with the image of such high density object(s) 
25 whilst still avoiding ihe generation of the above artifacts. 

Now therefore, according to one of its aspects the invmtion is characterized 
according to the characterizing part of Claim 1 . 

By itself US-A-4,590,558 discloses a method for removing objects from CT 
images, wherein an operator defines a **rub-out" region that encompasses the object to be 
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removed, whereafter the rub-out region is subjected to an averaging function. The operations 
by the operator clearly necessitate appreciable effort by a skilled worker who must carefully 
consider the possible location and shape of the high-density object, and on the basis thereof 
set the rub-out region. The prior art approach will nevertheless be prone to human and oth^ 

5 errors, and an automatic procedure would therefore be much preferred. 

Furthermore, US-A-6,094,467 requires, next to the standard imaging apparatus 
an additional hardware facility with many narrow-beam detection facilities for determining 
the extent, and in particular, the boundaries of bigfh attenuation objects and for thereby 
reducing the artifacts that the high attenuation objects would cause, without removing the 

10 high attenuation objects from the image. The additional plurahty of narrow beam detection 
facilities represent additional cost and additional control operations, and in consequence, the 
present inventors have imdertaken to derive all necessary information from a single 
measuring system. 

The invention also relates to a method being arranged for implementing the 
1 5 apparatus as claimed in Claim 1 , and to a computer program and to a computer program 
product compristQg instructions for controlling hardware for thereby being arranged for 
implementing the method as claimed in Claim 1. Further advantageous aspects of the 
invention are recited in dependent Claims. 

20 BRIEF DESCRIPTION OF THE DRAWING 

These and further aspects and advantages of the invention will be discussed 
more in detail hereinafter with reference to the disclosure of preferred embodiments, and in 
particular with reference to the appended Figures that show: 

Figure 1, a flow diagram of a procedure according to the present invention; 
25 Figures 2a-2e, the principle of suppression of the high-density objects; 

Figure 3, a 3D reconstruction image of metallic screws in a human knee 
envkoimient, executed whilst including the suppression and superimposing method for high- 
density objects according to the present invention. 

30 DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

Figure 1 illustrates a flow diagram of a procedure according to the present 
invention. In block 20, the procedure is started and the necessary hardware and software 
facilities are assigned. In block 22 the measuring apparatus executes a measuring scan. The 
prime considered technology is 3D-XRAY (3D-RX), but 3D-CT (computer tomography) 
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technology could be enhanced as well. In block 24, the system determines whether all 
intended scans have been executed. If no, the next scan is made after an appropriate rotation 
of the apparatus; if yes, the procedure proceeds to block 26. hi block 26, a ramp filter in the 
direction of rotation is s^plied to the projection. The setting of the ramp may be done once 
5 and for all, it may be done on the basis of statistical processing, such as based on tiie 

assumption titiat the area covered by the high density object(s) is generally small, or it may be 
done in a heuristic or even intuitive manner, or it could be effected by an operator person. 
Inasmuch as the high-density objects would generally cause a much greater attenuation that 
the tissues, tiie precise setting of the discrimination threshold is not critical. 

10 Next, in block 28, the filtered out high-density object(s) are segmented to get a 

3D reconstruction thereof. Next, in block 30 the high-density objects are suppressed m the 
original projection. Advantageously, they are then replaced by one or more gray values. This 
may be done by linear interpolation between the neighbouring pixels, by replacing each 
suppressed object by a single standard value, or by some other appropriate steps. In block 32, 

15 the ramp-filtered projections are then segmented a second time, but now without the high- 
density object(s). A straightforward approach would be to use exactly the same 
discrimination threshold as for the input values to block 28. 

Finally, in block 34, the segmented one or more higji-density objects from 
block 28 are superimposed on the segmentation result of the remainder. Finally, in block 36, 

20 the procedure is terminated, and the assigned facilities are relinquished again. The flow chart 
represents in various respects a simplification. For example, no escape procedure has been 
shown other than at the successfiil termination of the processing. Further, a trial and error 
procedure could be used for effectively setting the ramp threshold. 

Figures 2a-2e, illustrate the principle of suppression of the high density 

25 objects. The method will be described for 3D-RX. A 3D-imaging modality bases on a 
number of projections acquired diuing a rotational run with a motorized C-arm system. 
Figure 2a shows one simplified projection image. Therein, items 1 to 3 are high-density 
objects or so-called bullets used for registration of different data sets. After registration, the 
. bullets must be eliminated for a 3D-reconstruction in order to avoid artifacts. When 

30 measuring an intensity curve across the projection, at the position of the bullet a signal of 
high-intensity arises. This signal has to be suppressed and may be filled with gray values 
derived from the neighbourhood that comprises surrounding structures with lower densities. 

The suppression of the higji-density object(s) is effected by using ramp- 
filtered projections, which at the same time improves and simplifies the reconstruction. In 
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this respect. Figure 2a shows a projection image with bullets 1, 2, 3 clearly visible in the fonn 
of black circles. In a first processing step, all projections are ramp-filtered and used for a SD- 
reconstruction through back-projection with the well-known Feldkamp algorithm. Also, the 
bullets, or for that matter, other high-density objects, are segmented in the 3D-reconstraction, 
5 resulting in the picture of Figure 2b. The ramp-filtering results in sharper edges of the bullets, 
and the 3D-rec6nstniction enhances the contrast Through so acting, the segmentation of the 
high-density objects from the surrounding tissue or structures of their neighbourhood can be 
effected much better. Next, the bullets can be discrindbtxated by a simple threshold, so that in 
consequence only the bullets themselves are visible in the form of points with encircling 

10 rings. The latter represent the shadow of the bullet in question. 

Now, in a third processing step, the ranaining bullets are forward projected 
into ramp-filtered versions of the original projections in order to mark the search regions for 
the respective bullet borders, leading to the result illustrated m Figure 2c. After the above 
processing, the detected borders of the bullets will be forward projected into the original 

15 projections, resulting in the image shown in figure 2d. In the fourth processing step, the 
bullets are suppressed by substituting them by the gray values of the structures in their 
surroundmg neighbourhood, such as by linear interpolation between the respective entry 
point and the corresponding exit point of the bullet in question. Subsequently, a new 3D- 
reconstruction is performed through again applying flie Feldkamp algorithm with ramp- 

20 filtered projections. Thereby, the bullets are reduced or even eliminated, such also including 
removal of the artifacts caused by the bullets, leading to the image shown in Figure 2e. 

Figure 3 illustrates a 3D reconstruction image of metallic screws in a human 
knee environment, executed whilst including the suppression and superimposing method for 
high-density objects according to the present invention. The reconstruction clearly shows the 

25 knee joint made up of two major bones and the knee-cap, and also two screws used for 

clinically fixating bone parts to each other. Furthermore, clearly, no artifacts can be seen in 
the image. The result is attained through including the high-density suppression method in 
the 3D-reconstruction, The final image is reconstructed by matching the 3D-data set of the 
. bone structures with the 3D-reconstruction of the segmented screws. Furthermore, the gray 

30 values of the screws has been adapted to show both screws and bones at the same time. Such 
adapting may be represented by reducing the overall dynamic representation of the various 
objects, and in particular, of the high-density objects. Such simultaneous presentation allows 
a better diagnosis of the region between the screws and the surrounding bone structures- 
Artifacts caused by the high-density implants can be much reduced now. 
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The inventors have found that the results according to the present invention 
are much better when using ramp-filtered projections in combination with a 3D- 
reconstruction like througih the Feldkamp algorithm. The edges of the high density bodies 
will generally be much sharper and the contrast is improved. Thereby, a segmentation of the 
5 implants from the surrounding structures can be done much better. The segmentation can be 
easily performed with a simple threshold. This is an important aspect of the present 
invention, inasmuch as segmentation algorithnois are often qidte complex and thereby, time- 
consuming. 
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CLAIMS: 



1 . An apparatus for generating images of a human or animal body on the basis of 

3D-constructions from 3D-XRAY or 3D-Computer Tomography measurements, which 
bodies comprise both natural tissue and one or more high-density objects, said apparatus 
comprising a measuring faciUty for executing said measxirements, a distinguishing facility for 
5 distinguishing said one or more high-density objects and executing a separating procedure 
thereon for generating an improved image of regions of said natural tissue, 
said apparatus being characterized by: 

a ramp-filtering facility for applying a ramp filter in the direction of rotation to 
such various projection measurements and a back-projecting facility fed by said ramp- 
10 filtering facility for back-projecting the various so filtered projections into a 3D-volume 
reconstruction (Figure 2b); 

a segmCTLting facihty fed by said back-projecting faciUty for in said 3D- 
volume reconstmction segmenting said one or more high-density objects by a thresholding 
procedure and a forward projecting facility fed by said segmenting facility for executing a 
15 forward projection of the shadow(s) of the segmented one or more high-density objects onto 
the ramp-filtered projection (Figure 2c), whilst marking the borders of said one or more high 
density objects in the ramp-filtered back-projections; 

a suppressing faciUty fed by said forward projecting faciUty for suppressing 
said reconstructed one or more higji-density objects from the original projection 
20 measurementsand said suppressing faciUty is operative for executing an appropriate 

substitution of gray values derived from a physical neighbourhood of said one or more high- 
density objects instead of said one or more high-density objects in question. (Figure 2d); 

and a retro-coupling faciUty fed by said suppressing faciUty for executing a 
back-projection of the various filtered projections with corrected profiles through exclusion 
25 ' of said suppressed one or more high-density obj ects and outputting a reconstruction result 
OFigure 2e). 
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2. An apparatus as claimed in Claim 1, and furthermore comprising a 

superimposing facility fed by said forward projecting facility for receiving said one or more 
high-density objects for superimposing thereof onto said reconstruction result. 



3. An apparatus as claimed in Claim 1, and comprising adapting means for 

relatively adapting the gray values of said one or more high-density objects and said natural 
tissue in a predetermined gray value range to show both of them at the same time. 



4. A method for using an apparatus as claimed in Clsdm 1, for generating images 

of a human or animal body on the basis of 3D-constmctions firom 3D-XRAY or 3D- 
Computer Tomography measurements, which bodies comprise both natural tissue and one or 
more high-density objects, said method comprising the steps of executing said measurements, 
distinguishing said one or more high-density objects and executing a separating procedure 
thereon for generating an improved image of regions of said natural tissue, 
said method being characterized by comprising the steps of: 
applying a ramp filter in the direction of rotation to such various projection 
measurements and back-projecting the various filtered projections into a 3D-volume 
reconstruction (Figure 2b); 

in said 3D-volume reconstmction segmenting said one or more high-density 
objects by a Ihresholding procedure and executing a forward projection of the shadow(s) of 
the segmented one or more high-density objects onto the ramp-filtered projection (Figure 2c), 
thus marking the borders of said one or more high density objects in the ramp-filtered back- 
projections; 

suppressing said reconstructed one or more high-density objects firom the 
original projection measurements whilst executing an appropriate substitution of gray values 
derived firom a physical neighbourhood of said one or more high-density objects instead of 
said one or more high-density objects in question. (Figure 2d); 

and secondarily executing a back-projection of the various filtered projections 
, with corrected profiles and thereby without said suppressed one or more high-density objects 
(Figure 2e). 



5. A computer program comprising instructions for executing the method steps 

as claimed in Claim 4 through controlling an apparatus as claimed in Claim 1 . 
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6. A computer program product being embedded in a machine read-only tangible 

medium and containing instructions for executing the method steps as claimed in Claim 5 
through controlling an apparatus as claimed in Claim L 
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